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The chaotic ray dynamics inside dielectric cavities is described by the properties of an invariant 
chaotic saddle. The localization of the far-field emission in specific directions, recently observed in 
different experiments and wave simulations, is found to be a consequence of the filamentary pattern 
of the saddle's unstable manifold. For cavities with mixed phase space, the chaotic saddle is divided 
in hyperbolic and nonhyperbolic components, related, respectively, to the intermediate exponential 
(t < tc) and the asymptotic power-law {t > tc) decay of the energy inside the cavity. The alignment 
of the manifolds of the two components of the saddle explains why even if the energy concentration 
inside the cavity dramatically changes from t < tc to t > tc, the far- field emission changes only 
slightly. Simulations in the annular billiard confirm and illustrate the main results. 
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I. INTRODUCTION 

Motivated by technological applications, the emission 
from dielectric microcavities of different shapes has been 
the focus of detailed experimental investigations [U [2j |3l 
[10 and wave simulations [6, d H 10, 11, 11113 HI]. 
Comparisons between the measurements and the pre- 
dictions of the (chaotic) ray dynamics reveal an over- 
all good agreement, including detailed properties of 
the far- field emission ^1 These recent re- 

sults renew the interest on the classical (ray) dynam- 
ics in open chaotic systems per se, i.e., not only as the 
short-wavelength limit of the quantum (wave) descrip- 
tion [fs^lIlIIilIIHlIIlEolET]. 

The ray dynamics inside dielectric cavities is deter- 
mined by the laws of geometric optics: rays travel in 
straight lines between collisions at the boundary of the 
cavity, where they generically split in reflected and trans- 
mitted (refracted) rays with intensities given by Fres- 
nel's law. Far field emissions peaked in specific direc- 
tions have been surprisingly observed even in cavities 
where the reflected rays have (uniformly) chaotic dy- 
namics [Il|2l|3l[l|6l|8l[lll[Il|22]. Directionality 
in the far field and good confinement (high Q modes) 
are requirements for applications of microcavities as las- 
ing systems [1 . The following two recent results have 
proved to be crucial for a ray description of the far 
field emission: (i) Lee et al. [7 introduced the survival 
probability distribution of the intensities of rays inside 
the cavity. In strongly chaotic systems, this distribu- 
tion decays exponentially in time and numerical evidence 
was presented for a steady phase-space dependence of 
the probability distribution, independent of initial con- 
ditions (see Ref. [23^ for a detailed description, and also 
Refs. [7lll[l0l[l2l[19]). (ii) Schwefel et al [8 explained 
the observation of peaked far-field intensities by relating 
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the regions of high emission to the unstable manifold of 
periodic orbits close to the boundary of the region of to- 
tal internal reflection. This has been further investigated 
in Refs. [SliailOlinilTll. 

In this paper, the ray dynamics in dielectric cavi- 
ties is described using the ergodic theory of transient 
chaos 22 ESI EZ] • After properly taking into account 
the partial leak characteristic of dielectric cavities, the 
long-time properties of the chaotic dynamics are shown 
to be governed by an invariant set of the classical dy- 
namics, the so called chaotic saddle (CS), composed by 
all trajectories that never leave the cavity in both for- 
ward and backward times. The results from the theory 
of transient chaos are then applied, what gives the cor- 
rect theoretical framework for the description of chaotic 
optical cavities. In particular, a more general interpreta- 
tion of the results (i) and (ii) mentioned above becomes 
apparent: (i) the steady survival probability distribution 
introduced by Lee et al. in [7^ is equivalent to the con- 
ditionally invariant density pc [28l [29l |30] ; and (ii) the 
emission pattern is governed by the unstable manifold of 
the CS (not only of a single periodic orbit) along which 
the c measure concentrates [29j. The importance of the 
CS and its manifolds in quantum open systems has been 
recently recognized to explain the distribution of reso- 
nances [15] |T7l [19] and as the origin of a fractal Weyl's 
law [I6l [iSl [21] . Here, instead, I focus on the importance 
of the CS for the classical ray dynamics inside dielectric 
cavities, which are partially open systems. I find that 
the main physical observables (decay rate 7, emission 
pattern) can be obtained from properties of the CS. Fur- 
thermore, I argue how to extend these results to the case 
of generic cavities, where regions of regular and chaotic 
motion coexist in a mixed phase space. In particular, I 
show how a division of the CS in hyperbolic and nonhy- 
perbolic components [31, 32, 33 explains why even if the 
energy concentration inside the cavity changes in time, 
the far field emission retains its main properties. 

The paper is divided as follows. In Sec. [ll] the classi- 
cal ray dynamics and the standard description in terms 
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of the chaotic saddle is presented. The case of systems 



with mixed phase space is considered in Sec. |III| Sec. [IV 
presents numerical simulations on the annular billiard. 
Finally, the main conclusions are summarized in Sec. IVl 



II. CHAOTIC RAY DYNAMICS IN 
DIELECTRIC CAVITIES 

A. Classical ray dynamics 

Rays inside a dielectric cavity travel in straight lines 
between successive collisions at the cavity's boundary, 
where the ray generically splits in a reflected and a trans- 
mitted (or refracted) ray. The direction of propagation 
of the rays are determined by the angle with respect 
to the boundary's normal vector at the collision point. 
The reflected angle Or = is equal to the incident an- 
gle Oi, while the transmitted angle Ot is given by Snell's 
law as sin^T = nsin^/, where n is the ratio between 
the (constant) refractive indices inside and outside the 
cavity. The intensities of the rays after collision are 
given by Fresnel's law and total internal reflection oc- 
curs for p = sin 6 1 > sin 6c = l/n = Pc These are the 
well established laws of geometric optics. 

Assuming the validity of geometric optics, the dynam- 
ics of a ray is defined exclusively by its initial condition 
and the geometry of the cavity's boundary (parametrized 
by s). For simplicity, let us consider the case of two di- 
mensional cavities or billiards (the main results below 
remain valid for the three-dimensional case). The bound- 
ary's geometry defines a function M that maps one col- 
lision to the next M : (st^Pt) ^ (s^+i,pt+i). The map 
M preserves the area d/i = dsdp = dsd sin 0^ which es- 
tablishes the analogy to Hamiltonian systems. Below, 
the dynamics of maps M that have at least one chaotic 
component are considered. The discrete time t can be 
related to the actual time using the mean time between 
bounces ttA/Sc, where A is the area of the bilhard, S is 
the perimeter, and c is the speed of the ray. This is an 
approximation for individual rays [34 . 

The above description determines the dynamics in 
closed billiards. In order to introduce the escape through 
the transmitted rays (according to Snell's and Fresnel's 
laws), we consider that each ray has an intensity z, 
with it=o = 1. After each collision the intensity of the 
reflected ray it+i depends on i^, the angle Oj and on the 
polarization of the incident ray. For transverse magnetic 
(TM) and transverse electric (TE) polarizations, the re- 
flection coefficient R is given by the square of Fresnel's 
coefficients 



Rtm{0) = 
Rte{0) = 



sin{eT-Oi) 



sin{eT-\-Oi) 
tan(6'T+6'i) 



(1) 



for |sin(6>/)| < l/n = and R = 1 otherwise (to- 
tal internal reflection). The transmitted rays have an- 
gle Ot and intensity T = 1 — R. The region of the phase 



space —pc < p < Pc where T > 0, will be denoted as leak 
region I. The reinjection of transmitted rays into the 
cavity is neglected (it may occur in concave billiards). 

In summary, the full ray dynamics is given by 
(st.pt.it) ^ (s^+i,p^+i,zt+i), where M : (st.Pt) ^ 
{st-\-i,Pt-\-i) is an area preserving map defined by the ge- 
ometry of the billiard and the intensity z^+i = R{p)it 
decreases in time according to Fresnel's law Q. The en- 
ergy in one region O of the phase space at time t is given 
by the intensities it and density p{s^p^ t) of rays inside it: 

E{{s,p) eQ.,t)= / / i{s,p,t)p{s,p,t)dsdp. (2) 

The direction, position, and intensity of the rays emitted 
from the cavity can be computed by Fresnel's and Snell's 
law from E{s^p^t) inside /. 



B. Estimations based on the closed system 

Let us first consider the case of billiards where the 
dynamics of the closed map M is ergodic and strongly 
chaotic (e.g., uniformly hyperbolic) [26], leaving the 
generic case of systems with mixed phase space for 
Sees. Ill and IV For strongly chaotic systems, after a 
short transient time t* the fraction of rays that never 
entered / decays exponentially [26l [27] . The intensity i 



decreases with successive bounces inside /, and rays in / 
return typically exponentially fast to it [32 . Therefore, 
the total energy E{t) inside the cavity [i.e., considering 
Q in Eq. (l2|to be the full phase space] also decays expo- 
nentially \12\ [23] . This exponential decay is generically 
written as 



E{t) - (1 - rY = exp[ln(l - r)t], 



(3) 



where " ~ " indicates that both sides of the relation ap- 
proach a constant for long times, and the constant leak- 
age rate r corresponds to the transmitted energy per unit 
of time [26 . The escape rate 7 of Eq. ([3| is defined as 

7 = — ln(l — r) [~ r for small r]. (4) 

The ergodicity assumption for the closed system means 
that its phase space cannot be divided in two dynamically 
disjoint regions A and B with fi{A) > and p{B) > 0. 
Any initial density of rays po{s^p) = p{s^p^t = 0) con- 
verges (exponentially fast for strongly chaotic systems) 
to the natural density (constant in the phase space 
area dpds). When the leak / is small, a popular simpli- 
fying assumption is to consider the rays inside the open 
billiard at a given long time t to be distributed accord- 
ing to p^, i.e., according to the natural measure p of the 
closed billiard. Under this assumption, and taking into 
account that the transmission at time t occurs from in- 
side / according to T = 1 — i?, an approximation r* for 
the leakage rate r in Eq. Q can be calculated as 



T{e)dp= / / [i-R{e)]cosededs. 

Ji Jo J-Oc 



(5) 
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The leakage rate r was called the degree of leakage by 
Ryu et al. in Ref. [23 , where analytical expressions 
for approximation ([5|, using Rtm and Rte given by 
Eq. 0, 

were calculated. Ryu et al. show the in- 
teresting dependency r ~ l/v? that was verified nu- 
merically. Similarly, the ray dynamics described above 
has been successfully applied in cavities with different 
shapes [DlElllEllIllliaiTOllIllllIll]. It 
is interesting to compare these applications in optics 
to previous investigations involving other systems with 
leaks [281 |32l ESI |361 IMl [391 The main difference is 
the partial leak through Fresnel's law in the optical sys- 
tems (also present in acoustics [34]), in opposition to a 
complete escape assumed in the previous cases. However, 
as we will see below, once the intensity of the rays is prop- 
erly taken into account, a complete correspondence can 
be established. For instance, relation ([5| is a standard 
estimation [25l [26l [27] . Dependence on the position of / 
has been reported [SHl [38] and, as expected, Eq. ([5| is 
strictly valid only in the limit of r ^ 0, which is of little 
practical interest for optical systems since it corresponds 
to n ^ oo. However, it is important to note that ap- 
proximations based on closed system's properties, such 
as the one leading to Eq. ([5|, are often the only avail- 
able ones and lead to successful predictions [2 . In the 
next section, the analogy to leaked systems is deepened 
and fundamental results of the ergodic theory of tran- 
sient chaos (and chaotic scattering) are used to obtain a 
description of ray dynamics that fully incorporates the 
openness of optical cavities. 



ifold of the CS. A well defined escape rate 7 exists which 
is independent of the density of initial conditions po, pi'o- 
vided the support of po intersects the stable manifold of 
the CS. Relations between 7 and the properties of the CS 
(fractal dimension along the manifolds, Lyapunov expo- 
nent) have previously been derived [25} [26 } [27 ] . 

Let us see now how these results can be adapted to the 
case of optical cavities, where the leakage is only partial 
inside / and the dynamics involves not only the map M 
but also the decay of the intensity i. A natural definition 
of the CS is obtained replacing the condition of never 
escaping trajectories mentioned above by the condition 
that z = 1 for all times: 

{pes, scs) e CS 4=^ i{pcs, scs, t ±00) = 1. (6) 

In other words, the CS defined Eq. ^ is the same ob- 
tained considering a full leak and the standard definition 
of the saddle, i.e., replacing the partial leak in Eq. ([T|) by 
a Heaviside step function. It follows that: CS H / = be- 
cause escape takes place if (p, 5) G /. A definition similar 
to (|6| is not appropriate for the manifolds of the CS. For 
instance, initial condition inside / that converge to the 
CS for t ^ 00 and still have h^oo 7^ clearly deserve to 
belong to the saddle's stable manifold (similar argument 
for t ^ —00 holds for the unstable manifold). The stable 
(unstable) manifold of the CS is therefore defined by all 
points (s^p) CS for t +00 (—00), but attached to 
these points there is a manifold intensity i given by it^oo 
(it^-oo)- This suggests an alternative definition for the 
CS itself: 



C. Description in terms of invariant sets of the 
open system 

Transient chaotic motion is typical in systems with 
leaks and in naturally open systems showing chaotic scat- 
tering [24^ 25 , 26^ 27 . The escape of trajectories that 
remain a long time inside the system is governed by an 
invariant, nonattracting, chaotic set [24l|25[ [26l [27] . This 
set is composed by the trajectories that never leave the 
system, neither in forward nor in backward iterations of 
the map. The stable (unstable) manifold of this set is 
defined by all points that lead to this set in forward 
(backward) time. The term chaotic repeller is some- 
times used to denote this set [H [H [El [21 [25]. Be- 
cause billiards have a two dimensional phase space and 
are time reversible, having therefore stable and unstable 
manifolds, the term chaotic saddle (CS) is more appro- 
priate [271 [29]. For strongly chaotic ergodic systems, the 
CS has zero Lebesgue measure ii{CS) = (vanishing 
area of the phase space), and the support of the CS is a 
fractal set. The stable and unstable manifolds cross or- 
thogonally (angle bounded from zero) and are also of zero 
Lebesgue measure. Trajectories that survive for a long 
time inside the system necessarily have initial conditions 
close to the stable manifold of the CS, approach closely 
the CS, and leave the system through the unstable man- 



{Pcs, Scs) e CS ^ i{pcs, Scs, t ±00) > 0. (7) 

The CS defined using Eq. (ItI contains the CS defined 
using Eq. CS-^ C CS-^. Moreover, ah points in 
CS-Q but not in CS-(|6|, collide only a finite number of 
times inside / and necessarily belong to the intersection 
of the stable and unstable manifolds of CS-(|6|. There- 
fore, for long times, all rays with nonvanishing intensities 
inside the cavity will be governed by the CS (|6| and its 
manifolds. This justifies the choice of relation (|6| and 
shows that this CS also governs the energy decay from 
the billiard. 

With the above definitions of the invariant sets, let 
us characterize the escape from the cavity. The proper 
measure to describe this decay is the conditionally in- 
variant measure (c measure) djUc [28l [29l [30]. Intu- 
itively, the mathematically well defined c measure is ob- 
tained multiplying the survival density by a factor pro- 
portional to exp(7t) that compensates the decay of the 
Lebesgue measure fST. The conditionally invariant den- 
sity pc concentrates along the unstable manifold of the 
CS [2^ and is the only attractor for typical initial densi- 
ties po [28l [30] . An important property of the c measure 
is that it converges to the natural measure for small leak 
regions p{I) [W, W, which justifies the approxima- 
tion used to obtain Eq. ([5|. However, typical dielectric 
cavities have n < 10 and the approximation of small leak 
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is violated. This means that the dynamical properties 
derived from the geometry of the (closed) billiard are 
not a good approximation to the dynamics of the optical 
(open) billiard. In analogy to the calculation in Eq. ([s]), 
but this time using the precise distribution inside the 
cavity given by the leakage rate r in Eq. ([3| is given 
by the c measure of the leak / [2811321 [37j 

r /"l r-\-Oc 

r= T{0)d/j.c = [1 - R{0)]pc{0,s) cos OdOds. 

Ji Jo J-e^ 

The escape rate 7 is obtained from Eq. Q. For small / 
(large n) pc approaches and r in (|8| approaches r*, 
obtained in ([5). Typically r* overestimates r [33 . 

It is remarkable that many of the results of the theory 
of transient chaos presented above have been recently 
rediscovered in the analysis of optical cavities. For in- 
stance, the steady state of the survival probability dis- 
tribution, introduced by Lee et al. [7 and mentioned as 
point (i) in the Introduction, is equivalent to the condi- 
tionally invariant density pc introduced by Pianigiani and 
Yorke thirty years ago [28] (see Refs. [29l |30] for a recent 
review). The independence of initial ensembles reported 
in Ref. [7] is related to the existence of the invariant CS 
or, equivalently, to the fact that pc is the only attr ac- 
tor for typical po's. The density pc concentrates along 
the unstable manifold of the CS [29], which presents the 
characteristic filamentary pattern inside /. This explains 
the peaked distribution of the transmitted rays responsi- 
ble for the emission from the cavity. When a short time 
periodic orbit exists close to the critical line p = Pc^ the 
unstable manifold of this orbit (which also belongs to 
the CS) will be parallel to the manifold of the remaining 
part of the CS because both manifolds do not intersect. 
This corresponds to the observation by Schwefel et al. [8] 
described as point (ii) in the Introduction. 



III. SYSTEMS WITH MIXED PHASE SPACE 

The results described so far can be rigorously applied 
only for a limited class of strongly chaotic systems. It can 
be argued that it is a technical problem to extend these 
demonstrations to a larger class of nonuniformly hyper- 
bolic systems, where no deviation of the exponential de- 
cay have been numerically detected. However, cavities 
with generic boundaries typically have a mixed phase 
space: coexisting with regions of chaotic motion there 
are regions of regular motion, e.g., Kolmogorov-Arnold- 
Moser (KAM) islands. Around these regions there is a 
sticky region whereto chaotic trajectories get partially 
trapped, introducing a power-law like decay of the sur- 
vival probability. Rigorously, the results of the theory of 
transient chaos mentioned above either do not apply or 
become trivial. This type of cavities have been consid- 
ered in Refs. |3| [61 |8j and without further justifications 
it is not clear how the previous results can be extended 
to this case. 



In this section, I show how, despite the mathematical 
difficulties, in practice the formalism of the CS and its 
manifolds do apply to billiards presenting a divided phase 
space. The basic observation is that typically a well de- 
fined exponential decay of the survival probability exists 
for intermediate times, where in practice the previous re- 
sults can be applied [31, 32 . Below, billiards containing 
a large chaotic component and no KAM islands in the 
border of / are considered. In this case, similar to the 
decay of trajectories, the decay of energy ([2| from the 
chaotic component shows a transition from exponential 
to asymptotic power-law [23l l32l [35l [39 l 140 } [IT] 



Ae 
Ae 



-It 
-It 



-Bt-' 



for t > r, 

for t > ta^ 



(9) 



where is proportional to the inverse of the negative 
Lyapunov exponent of the saddle, t^ > t* is the short- 
est time rays in the sticky region need to reach /, and 
A^-itc^ ^ Bt~^. The physically relevant time is the 
crossover time tc > t^ between the exponential and 
power-law decays in Eq. It is defined as Ae"^^^ = 
Bt~^ = E{tc)/2, i.e., for t > tc the power-law decay dom- 
inates. In Ref. [32 it was shown that tc ~ I/7 and sug- 
gested that exponential and power-law regimes in Eq. (|9| 
can be related, respectively, to hyperbolic and nonhyper- 
bolic components of the CS, as first suggested in Ref. [31j. 
The nonhyperbolic component consists of the border of 
KAM islands (and by other marginal stable orbits) while 
the hyperbolic part is away from the sticky regions and 
resembles the CS described in Sec. Ill Cl Initial conditions 
uniformly distributed touching the KAM islands reduces 
the exponent a in Eq. ([9| by 1 (see Ref. [41 and refer- 
ences therein). Considering this effect, and because ([9| is 
a survival probability distribution, typically < a < 1. 
The exact (finite time) value depends on the properties 
of the nonhyperbolic region (KAM island). 

In terms of the energy inside the cavity, Eq. ([9| indi- 
cates that while for t < tc energy concentrates strongly 
in the hyperbolic component of the CS and its mani- 
folds, for t > tc the energy concentrates in the nonhy- 
perbolic component of the CS close to the KAM islands. 
In principle, one could expect also a dramatic change in 
the emission pattern from t < tc to t > tc. However, 
as emphasized in Ref. [32], one important difference be- 
tween the hyperbolic and nonhyperbolic components of 
the CS is that their manifolds attract and repel expo- 
nentially and sub-exponentially respectively. Therefore, 
rays slowly approach and slowly escape the nonhyper- 
bolic component of the CS through the hyperbolic compo- 
nent. In particular, when KAM islands are away from /, 
rays that approached KAM islands will come close to 
the stable/unstable manifolds of the hyperbolic compo- 
nent of the CS before escaping. More precisely, the un- 
stable manifold of the nonhyperbolic component of the 
CS aligns to and follows closely the unstable manifold 
of the hyperbolic component. From this qualitative de- 
scription we expect that, even if the energy concentrates 
for long times t > tc on the nonhyperbolic component of 



FIG. 1: (Color online) (a) Annular billiard defined by two eccentric circles, with control parameters q — 0.65 and 6 — 0.3. 
Lines indicate one stable periodic orbit (at = 0) and one MUPO (7,1) (polygonal form) [42^. Trajectories that do not cross 
the dotted circumference of radius q-\- S belong to the whispering gallery, (b) Phase space representation of the closed billiard 
shown in (a), obtained using a Poincare surface of section at the outer boundary (parametrized by s = G [0, 27r]). The stable 
periodic orbit in (a) is at the center of a chain of KAM islands. Black dots mark 8000 iterations of a single chaotic trajectory. 
The horizontal dashed lines denote pc = ±1/^ for n = 4, as used in the following figures. 



the CS, the decay towards / will not dramatically change 
from t < tc to t > tc. This prediction is partially con- 
firmed in the next Section, where numerical simulations 
of a specific system are presented. 



IV. NUMERICAL SIMULATIONS 

A. The annular billiard 

Numerical simulations of the ray dynamics are per- 
formed in the annular billiard, composed by two circles 
(of radii q < 1 and R = 1) placed eccentrically (distance 
5 between centers), as depicted in Fig. [l^. The closed 
phase space shown in Fig. shows a large chaotic sea 
coexisting with three regions of regular dynamics: the 
chain of KAM islands and two whispering gallery regions 
close to the outer boundary. Whispering gallery regions 
are typical for billiards with concave shape [4,|6l[T4]. In 
the case of the annular billiard an infinite number of fam- 
ilies of marginally unstable periodic orbits (MUPOs) ac- 
cumulate outside of the whispering gallery [42 , whereto 
the chaotic trajectories stick with a = 1 in Eq. ^ [41^ . 



B. Time dependence 

Consider now that the larger circle in the annular 
billiard (Fig. [l^) has refractive index n, the outside 



space n = 1, and the boundary of the inner circle is a per- 
fect mirror. In this case, emission according to Eq. ([T]) is 
possible through the outer border when \p\ < \pc\ = 1/n, 
i.e., the leak region / corresponds to the horizontal stripe 
between —pc <p < Pc in the center of the phase space 
shown in Fig. [l[b). Typically, 10^ TM polarized rays 
(similar results are expected for TE polarization) are 
started uniformly distributed according to the Lebesgue 
measure (area of the phase space) inside the chaotic com- 
ponent. This means that no rays are started inside the 
KAM island or whispering gallery. The temporal decay 
of the total energy is presented in Fig. [2] for different val- 
ues of the refractive index n. It confirms the existence 
of an intermediate exponential decay and an asymptotic 
power-law decay, as described by Eq. ([9|. Different ex- 
ponents can be identified for n = 2 and 4: 1 < a < 2 
for intermediate times, related to rays that approach the 
KAM island (started away from it), and < a < 1 re- 
lated to rays started already close to the KAM islands. 
For n = 2, 4, and 10 the energy decay converges precisely 
to the same curve for t oo, because the nonhyper- 
bolic component of the CS is the same in all these cases. 
For n = 1.5 this is not true because the KAM island 
chain is (partially) inside /. Hereafter the representative 
case n = 4 is chosen, for which the generic exponen- 
tial and power-law decays are clearly visible in Fig. [l] 
and tc = 80. 

Let us see now how the energy is distributed inside the 
billiard at a given time t. Figure |3] shows the distribution 
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FIG. 2: (Color online) Fraction of initial energy inside the bil- 
liard shown in Fig. [l] as a function of time. Different curves 
correspond to the refractive indices n = {10,4,2,1.5} (from 
top to bottom), (a) Linear-log and (b) log-log scale. The 
dashed lines correspond to an exponential fitting for short 
times and tc indicates the transition time to a power-law. 
Rays were started uniformly distributed in the chaotic com- 
ponent. 



projected to the p axis for different times t. The total 
energy in all cases is normalized and the rays are started 
uniformly in the chaotic component. For a fully chaotic 
system this distribution converges to the conditionally 
invariant density pc (see Sec. Ill) projected on the p- 
axis. Instead, Fig. |3] shows a concentration of the energy 
close to the KAM island and whispering gallery. This is 
a consequence of the power-law escape of rays in these 
sticky regions, in opposition to the exponential escape for 
rays away from these regions. Apart from a rescale due 
to this decay. Fig. [3] suggests that the distribution away 
from islands also follows a similar pattern (see, e.g., the 
vertical arrows in Fig. [3| for all times t < tc and t > 
tc. This is in agreement with the interpretation given 
at the end of Sec. |III| and will be further investigated in 
Sec. lTVDl 
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FIG. 3: (Color online) Energy along p at time t normalized 
by total energy at time t, for t = 0, 1, 40, 80, 160, and 300 (top 
to bottom). Refractive index n = 4 leads to pc — 0.25 and 
= 80 (see Fig. [2}. Rays were started uniformly distributed 
in the chaotic component [outside the KAM island). Vertical 
arrows indicate minima of the distribution that systematically 
appear for all times. The dotted line indicates the reflection 
coefficient Rtm in Eq. ([T]), multiplied by 10 for clarity. 



C. Chaotic saddle 

In this section we will see how the division of the 
chaotic saddle in hyperbolic and nonhyperbolic compo- 
nents proposed in Sec. |III| apply for the annular billiard 
considered here. The nonhyperbolic component of the CS 
is composed by the border of the KAM islands, the bor- 
der of the whispering gallery, and the families of MUPOs. 
Figure [3] shows that the energy concentrates in this region 
for long times. Regarding the hyperbolic component, a 
simple and efficient method to obtain a visualization of 
this zero measure set was proposed in Ref. [27] (p. 201). 
It is based on the observation that (most) trajectories 
that survive until some long time t"^ were close to the sta- 
ble manifold of the CS at time t = 0, were close to the CS 
at time t = t^ and were close to the unstable manifold 
at time t = t^ . This method is expected to work for the 
hyperbolic component of the CS if <C ff" < tc and if ini- 
tial conditions are selected away from sticky regions, i.e., 
non-uniformly in the chaotic sea. The results achieved 
for the annular billiard are shown in Figs. [4]^b) and[4]^c). 
Comparing to the closed system's phase space shown in 
Fig. Qa), it is evident that non-trivial structures were 
created by the openness of the system. Furthermore, in 
opposite to the density of the closed system [support 
shown in Fig.jlja)], the conditionally invariant density pc 
of the open system [support shown in Fig. ^c)] is not 
smooth. 



D. Energy distribution and emission 

Finally, the energy distribution inside the cavity and 
the far field emission are investigated. Considering the 
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FIG. 4: (Color online) (a) Closed system's phase space [mag- 
nification of Fig.[T]^b)]. (b) Hyperbolic component of the CS, 
and (c) its unstable manifold for n = 4. The method of 
Ref. 123 with = 30 < tc = 80 was employed, with rays 
started away from the regions of regular motion. Notice that 
the unstable manifold enters / in (c), a sign of rays leaving 
the cavity (only points with i > 0.01 are plotted). Due to the 
time reversible symmetry and the spatial symmetry of the an- 
nular billiard, the stable manifold of the CS can be obtained 
from the unstable one by (s,p) i-^ (27r — s,p). 



time dependency obtained in Sec. |IVB[ where tc = 80 
was found for n = 4, two times are considered: t = tc/^ = 
40, for which the decay is exponential and the hyperbolic 
component of the CS dominates, and t = 2tc = 160, for 
which the decay is algebraic and the nonhyperbolic com- 
ponent of the CS dominates. The relative (apart from 
the overall decay) phase-space distribution of energy at 



these times are expected to be representative for all times 
in the exponential and power-law decays (e.g., the nor- 
malized distribution for t ^ oo is expected to resemble 
the one at t = 160). The numerical results are shown in 
Fig. m The two fi gures in the upper row confirm that 
the energy shifts from the hyperbolic component of the 
CS and its unstable manifold, to the nonhyperbolic com- 
ponent of the saddle. The two figures in the bottom row 
are magnifications of the upper figures in the leak region, 
which is the relevant region for emission purposes. The 
differences between t = 40 and t = 160 are much less dra- 
matic in this case. Comparing to Fig. Qc), we see that 
the energy is non-vanishing along the unstable manifold 
of the (hyperbolic component of) the CS, and only the 
relative intensity is (slightly) changed in time. These dif- 
ferent intensities are enough to change the far-field emis- 
sion as shown in Fig. |6] An overall similar pattern is 
observed for both times, but with different intensities at 
different emission angles. These results are in agreement 
with the interpretation at the end of Sec. |III| that the 
unstable manifolds of the hyperbolic and nonhyperbolic 
components of the CS are close to each other inside /. 
However, the numerical simulations presented here sug- 
gests that differences in the intensities inside / lead to far 
field emissions that are similar but not identical for t < tc 
and t > tc. 



V. CONCLUSIONS 

This paper describes the chaotic ray dynamics in di- 
electric cavities in terms of the theory of transient chaos. 
After properly taking into account the partial leakage 
introduced by Fresnel's law, the long time escape and 
emission are found to be governed by a chaotic saddle 
(CS) [21 EH EH [27]. The energy inside the cavity is dis- 
tributed according to the c measure [28 , that is non-zero 
along the unstable manifold of the CS [29 . In the generic 
case of cavities showing mixed phase space, it is useful to 
consider a division of the CS in hyperbolic and nonhy- 
perbolic components [3TJ [32] , related to the exponential 
{t < tc) and power-law {t > tc) decays of energy inside 
the cavity. For longer times {t > tc) the energy concen- 
trates in the nonhyperbolic component of the CS, but 
the emission shows only small differences from shorter 
times {t < tc). This is explained by the fact that, when 
the nonhyperbolic regions are away from the critical line, 
the unstable manifolds of the hyperbolic and nonhyper- 
bolic components are close to each other inside the leak 
region. 

From the point of view of previous investigations in 
dielectric cavities, these results provides an unifying and 
more general description of experiments and simulations. 
For instance, the localization in the far field emission 
is shown here to be related to the unstable manifold of 
the CS and not of a single periodic orbit as previously 
reported, what is expected to be of practical relevance in 
cases where no simple periodic orbit exists close to the 
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FIG. 5: (Color online) Energy density in the phase space of the annular billiard for n — 4 (jpc — 0.25) at two different times: 
Left t — 40 = tc/2; Right t — 160 = 2tc. Bottom row shows magnifications of the top row in the leak region. Rays were started 
uniformly distributed in the chaotic region. 
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FIG. 6: (Color online) Far field intensity transmitted in the 
direction p — (f) — Ot measured from the center of the annular 
billiard at times t = tc/2 = 40 (dashed line) and t = 2tc = 160 
(solid line). Rays were started uniformly distributed in the 
chaotic region and n = 4. 



critical line. 

From the point of view of dynamical systems theory, 
the novel aspects of this paper are twofold: (I) it takes 
the partial leak of optical systems explicitly into account 
in the definition (of the unstable manifold) of the OS 
(Sec. II C); and (II) it shows how systems with mixed 
pha se space can be effectively described (Sees. |III| and 
[lV[ ). The success of point (II) in the description of phys- 
ically relevant quantities, such as the intensity of the far 
field emission, suggests that future works in nonhyper- 
bolic Hamiltonian systems should consider in more detail 
and more rigorously the division of the chaotic saddle in 
hyperbolic and nonhyperbolic components [31] |32] . Alto- 
gether, it is also interesting to see that many of the funda- 
mental concepts of transient chaos theory (e.g., the OS, 
its unstable manifolds, and c measure) achieve an exper- 
imental concreteness in optical systems. This is perhaps 
only comparable to the case of fluid dynamics, where the 
unstable manifold of the CS can be visually observed in 
scattering systems, specially when a constant injection or 
activity of tracers compensate the natural decay through 
the fluid flow EH |38l Hal . 



Similar to the case of fluids, experiments and applica- 
tions in dielectric microcavities are rarely a simple de- 
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cay of rays that have been excited at some time t = 0, 
as considered in the model of this paper. Instead, en- 
ergy is constantly pumped to the system from outside 
through some gain in the medium. The relevance of the 
rays that survive for long time inside the cavity by ap- 
proaching the CS is that their intensities are enhanced 
through the input of energy. Accordingly, in the wave 
picture, the (high-Q) lasing modes concentrate in the re- 
gion of the CS [15j. The details of these experimental 
mechanisms, and how they could be included in the ray 
picture, are beyond the scope of this paper. However, 
generally one can think that these different mechanisms 
translate in a minimum confinement time t^, and only 
rays confined for times t > to are relevant. In this pic- 
ture, the results of this paper suggest that depending 
whether to < tc or to > tc the energy inside the chaotic 



component of the cavity changes dramatically: it will be 
concentrated mainly in the hyperbolic (to < tc) or non- 
hyperbolic {to > tc) component of the saddle (in optical 
microcavities to ^ tc). However, the emission pattern is 
similar in these two cases because it is determined by the 
unstable manifold of the CS inside /. A similar emission 
pattern is also expected when the relevant modes are in- 
side the KAM islands and tunnel the (dynamical) barrier 
that separates them from the chaotic sea. 
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